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^ ' Abstract 

iu: 

^ ■ Recent experiments by the HERMES group at HERA are measuring semi- 



inclusive electroproduction of pions from deuterium. We point out that by 
comparing the production of 7r+ and vr^ from an isoscalar target, it is possi- 
ble, in principle, to measure charge symmetry violation in the valence quark 
distributions of the nucleons. It is also possible in the same experiments to 
obtain an independent measurement of the quark fragmentation functions. 
We review the information which can be deduced from such experiments and 
show the "signature" for charge symmetry violation in such experiments. Fi- 
nally, we predict the magnitude of the charge symmetry violation, from both 
the valence quark distributions and the pion fragmentation function, which 
might be expected in these experiments. 
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I. INTRODUCTION 



It has long been been recognized that charge symmetry or isospin symmetry is respected 
in the strong interaction to a high degree of precision. Charge symmetry violation [CSV] in 
nuclear physics is generally good to within about 1% of the strong amplitude [0,0. Therefore, 
in most analyses of strong interactions, it is reasonable to assume the validity of charge 
symmetry. 

Charge symmetry has always been assumed in defining quark distribution functions. It 
is so common in quark/parton phenomenology that it is frequently not even mentioned as 
an assumption. Charge symmetry reduces by a factor of two the number of independent 
quark distributions one must define, and until recently there has been no compelling reason 
to suggest charge symmetry violation. On the other hand, there were no precise tests of 
charge symmetry in parton distributions, either. 

At the nucleon level, charge symmetry appears as the symmetry under interchange of 
protons with neutrons. At the quark level, charge symmetry involves interchange of up and 
down quarks. Recently, the question of charge symmetry at the quark level has become of 
great interest. At the present time, we have a good quantitative knowledge of the valence 
and sea quark distributions in the nucleon. Furthermore, it is now possible to perform 
experiments whose precision, and range in x, is sufficient to probe "small" components 
of quark distributions, and hence to test quantities like flavor symmetry of nucleon sea 
distributions, and charge symmetry of parton distributions, which were not accessible to 
direct experimental test in the past. 

Interest in this general field has been sparked by the experimental evidence of violation 
of the Gottfried sum rule reported by the New Muon Collaboration (NMC) @]. Although 
this has widely been quoted as evidence of SU(2) fiavor symmetry violation [FSV] in the 
proton sea 0- (i.e. u^{x) ^ d^{x)), an alternative explanation of that result could be 
from charge symmetry breaking in the proton sea, as has been pointed out in Refs. (it 
could also arise from a linear combination of CSV and FSV effects in the nucleon sea). 
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Initial theoretical investigations of charge symmetry in valence quark distributions have 



been carried out by Sather |T0|, and by Rodionov et al. |TT]]. Following the approach of Ref. 

, they calculated valence quark distributions for the proton and neutron in a quark model 
and in the MIT bag model, respectively. In both models the "majority" quark distributions 
(i.e., 'U^(x) and satisfied charge symmetry to within about 1%, while the "minority" 

quark distributions {dP{x) and were predicted to violate charge symmetry by 5% or 

more at large x. More recently, Londergan et al. |]13| proposed that pion-induced Drell-Yan 
processes could be used to probe these CSV effects in nucleon parton distributions. 

In this note, we point out that semi-inclusive pion production, from lepton deep inelastic 
scattering on nuclear targets, could also be a sensitive probe of CSV effects in nucleon 
valence distributions. Several authors have previously discussed the possibility of using semi- 
inclusive scattering processes to address various physics issues [p!^ -[16[|. In fact, the HERMES 



collaboration at HERA [|T^ is currently taking experimental data on semi-inclusive pion 
production from hydrogen and deuterium. 

The organization of our paper is as follows. In Sec. II, we review the formalism for semi- 
inclusive hadron production at high energies. We show that for an isoscalar nuclear target, 
one can test charge symmetry violation in both the valence quark distribution functions, and 
in the fragmentation function. To first order in these small effects, there is an experimental 
combination for which these two quantities separate. In Sec. Ill, we derive predictions for 
ratios of vr"*" and vr" semi-inclusive electroproduction from deuterium. We show both the 
qualitative and quantitative "signatures" for charge symmetry violation in these reactions. 
We present our conclusions and suggestions for future experimental work in Sec. IV. 

II. CHARGE SYMMETRY TESTS IN SEMI-INCLUSIVE HADRON 

ELECTROPRODUCTION 
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A. General Formulas for Semi-Inclusive Reactions 



In the quark/parton model, the semi-inclusive production of hadrons in deep inelastic 

lepton scattering from a nucleon is given by 

1 da%{x,z) _ j:ieUf{x)D^{z) ^ 
aN{x) dz J2iefq^{x) 

In Eq. (|l|), qf{x) is the distribution function for quarks of flavor i, and charge Cj, in the 

hadron as a function of Bjorken x. D^{z) is the fragmentation function for a quark 

of flavor i into hadron h. The fragmentation function depends on the quark longitudinal 

momentum fraction z = Eh/u, where and u are the energy of the hadron and the virtual 

photon respectively. We write the numerator in Eq. (|lD as A^^'^ = X^i gf^(a;)-Df (-z), so the 

quantity A^^'' represents the yield of hadron h per scattering from nucleon A^, as a function 

of z and x. In terms of these quantities, semi-inclusive production of a charged hadron from 

a proton can be described by 

m^{x,z) = ^u^{x)Di{z) + ^u^{x)D^{z) + ldP{x)D^iz) 

+ ld^{x)Dl{z) + ls^{x)Df{z) + ^s^{x)Df{z), (2) 

where Df{z) are the fragmentation functions for a quark (or antiquark) of flavor i into 
positively or negatively charged hadrons. 

Charge conjugation invariance implies that = and D"^ = Dj. By making 
the additional assumption of charge symmetry, the fragmentation functions for pions from 
quarks will obey the relations 

Dnz)=Df{z) 

Dl\z)=Dl{z) (3) 

Later, we will consider the possibility of charge symmetry violation in the fragmentation 
functions. 

We want to find an experimental signature for charge symmetry violation in valence 
quark distributions of the nucleon. We will therefore derive expressions for 7r+ and vr" 
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electropro duct ion on an isoscalar nucleus. In this paper our expressions specifically refer to 
a deuteron target, although our results can be extended to any isoscalar target. 



B. Experimental Extraction of Fragmentation Functions 

Levelt et al. \TA\ derived a useful expression by which the ratio of fragmentation functions 



can be extracted from leptoproduction of charged pions on protons and neutrons. They 

proposed measuring the quantity 

_ [NP""* (x, z) - A^"'^^ (x, z) + NP^~ (x, z) - AT"^" (x, z)] 
^' ~ [NP-^* (x, z) - N''-^^ (x, z) - NP""' (x, z) + A^"'^" (x, z)] ^ ^ 

In the quark/parton model, R{x, z) has the form 

R(x z) - ^ ("iM^lM^ (DM + DM) 
^ ' ^"5 [r{x) + fix)) Wiz)-D~iz))^ 

with the definitions 

r(x) = uP{x) - (F{x) 

f{x) = -uP{x) + dP{x). (6) 

Integrating both the numerator and denominator of -R(x, z) over x gave the result 

Q{z) = dx [A^P"^ (x, z) - A^""^ (x, z) + NP^~ (x, z) - A^""" (x, z)]/ 
Jo 

( / rfx[A^P"^(x,z) - A^""''(x,z) - A^^'""(x,z) + A^"""(x,z)]) 



^9° ( D+iz)+D-iz) \ 

5 ^ W{z)-D-{z)) ^ ^ 

In Eq. d^), Sg is the experimental value for the Gottfried Sum Rule 0. 

The quantity Q{z) is a function of the ratio of the "favored" and "unfavored" fragmen- 
tation functions, -D^ (2;) and D~{z), respectively. They are so named because production of 
a positively (negatively) charged hadron will preferentially occur from the up (down) quark 
in the nucleon. The fragmentation functions have been extracted by the EMC group |T8|,|19[, 
but an independent measurement of this can be obtained in pion leptoproduction. Feynman 
and Field pOl suggested a form for the ratio 



D-(z) 



1 - z 



Dt{z) l + z 

If Eq. (I) is correct, then measurement of Eq. (0) would give 
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(8) 



Q{z) 



5z 



S, 



G 



(9) 



The quantity Q{z) of Eq. (0) requires charged hadron production from both protons 
and neutrons. For charged pion leptoproduction on an isoscalar target (e.g., deuterium), a 
useful quantity is A^^'^ — N^'^^ . For an isoscalar target, the contributions from the sea will 
exactly cancel, and the result will depend only upon the valence quark contributions. From 
Eq. (^, we can show that this quantity can be written 



N 



{x,z)-N'''' {x,z) 



^ Kix) + Kix)] - i K(x) + <(x)]) [dUz) 



D-(z] 



(10) 



In Eq. ([T0|), we assume the validity of the impulse approximation, i.e., 



X, z 



(x, z). Note that the valence quark distributions appearing in Eq. (|10D are 
the quark distributions in the deuteron and not the free nucleon valence quark distributions. 
Integrating this quantity over all x, we obtain 



dx 



N^^^ix^z) 



N 



D-n 



Dl\z)-Dl{z) 



(11) 



The quantity Nd{z) is proportional to the difference between the "favored" and "un- 
favored" fragmentation functions, D^^ (z) and (z), respectively. If we assume the 
Feynman-Field parameterization for the ratio of fragmentation functions, then measure- 
ment of Eq. ([TT]) would give 

2z 



Nd{z) 



l + z 



Dl\z) 



(12) 



This could then be compared with the fragmentation functions extracted from the EMC 
group, or by the method suggested by Levelt et al. which requires comparison of pion 
leptoproduction on both protons and neutrons separately. 
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The fragmentation functions have previously been measured in deep inelastic muon scat- 
tering by the EMC collaboration [p!8| , p!9| . In Fig. [la] we show the favored and unfavored 
fragmentation functions for quarks into charged pions, as measured by EMC. In Fig. [lb] 
we plot the experimental ratio of fragmentation functions A{z), and compare it with the 
Feynman-Field parameterization of Eq. (^). For moderate values of z, the Feynman-Field 
parameterization is rather accurate. For the largest values of z, the experimental errors are 
large, but the experimental results appear to be systematically larger than the Feynman- 
Field predictions. 



C. Charge Symmetry Violation in Semi-Inclusive Pion Production 



We want to measure charge symmetry violation in the valence quark distributions. For 
pion electroproduction on an isoscalar target, (such as the deuteron) this can be achieved 
by noting that (assuming charge symmetry) the 'favored' production of charged pions from 
valence quarks are related by 

iVf,f (x,z)=4iVf,f (x,z). (13) 

That is, for 7r+ (tt^) production, the "favored" mode of charged pion production is from the 
target up (down) quarks. Since the semi-inclusive reactions are proportional to the square 
of the quark charge, there is a relative weighting of 4 for vr"^ production. 
We therefore propose measuring the quantity R^{x,z), defined by 



R^ix.z) 



N^''+{x,z) -ND^-{x,z) 
5A{z) {A + A{z))6D{z) 



+ 



1 - A{z) 

l + A{z) 
1 - Aiz) 

AJz) 



3[l-AizW 
4 {5d{x) - Su{x)) + 15 {uP{x) + (F(x)) 



3 [u^{x) + d^x)] 



[s{x) + s{x)] 



^ 1 - A{z) [u^{x) + d^{x)] ^^^^ 
The quantity contains the charge symmetry violating quark distribution functions 
defined by 



Sd{x) =dP{x) -u'^ix) 
6u{x) = vF{x) - 



the strange/favored ratio of quark fragmentation functions 



and the charge symmetry breaking fragmentation functions, 

SD{z) = . (15) 

We can more cleanly separate the x and z dependence by multiplying by a z- 
dependent factor, e.g. 

R^{x,z)^^^^^^R{x,z) 

5A{z) (4 + A{z))6D{z) ^ 4 [6d{x) - 6u{x)] 



+ 



1 + A{z) 3(1-A2(2)) ' 3[ul{x) + d^.{x)] 
5 {uP{x) + dP{x)) + As{z) [s{x) + s{x)] /(I + A{z)) 
[uv{x) + dv{x)] 

Rfiz) + R^svi^) + RLix,^) , (16) 



where 



R?(z 



5A{z) {A + A{z))5D{z) 



' 1 + A(z) 3(1- A2(z)) 



4 [5d{x) - 5u{x)\ 
3 {ul{x) + (R{x)y 



5 {u^{x) + d^jx)) + A,{z) [s{x) + -s{x)] /(I + A(2:)) 
^'^'^^^'"'^ [m?(x) + d?(a;)] ■ ^ 

Eq. ([1^) is obtained from Eq. ( p^ by multiplying by the experimentally measured fragmen- 
tation functions. In Eqs. (|1^, (|16|) and (p!7D, we have expanded to first order in "small" 
quantities. These are: the CSV nucleon terms, 5d{x) and 5u{x)\ the CSV part of the frag- 
mentation function, 5D{z)\ and the sea quark distributions (Eq. ( [T6| ) is only valid at large x 
where the ratio of sea/valence quark distributions is small). We have neglected the CSV part 
of the "unfavored" fragmentation function. In the region of interest (moderately large z) the 



unfavored fragmentation function will be considerably smaller than the favored term, and 
consequently the unfavored CSV term should be proportionately smaller than the favored 
CSV term. 

The quantity R^{x,z) separates into three pieces. The first piece, Rf{z), depends only 
on z, as is shown in Eq. (|T^. It contains a small part which is proportional to the CSV part 
of the fragmentation function. The dominant piece of Rj' (z) has the form 

5A{z) ^, 5(l-z) 

TTW) 2^' ^ ^ 

where the relation in Eq. ([T^) follows if we adopt the Feynman-Field parameterization. For 
most values of z, this term should decrease monotonically as z increases (although, at the 
largest values of z, it may increase, as suggested in Fig. [lb]). The second term, R^gy{x) 
depends only on x, and is proportional to the nucleon CSV fraction (relative to the valence 
quark distributions). The term R^^^{x,z) is proportional to the sea quark contributions. 

Experimentally, one needs to measure accurately the x-dependence of R^{x, z) for fixed 
2;; in this case the z-dependent term will be large (of order one) and constant. The sea quark 
contribution will be large at small x, but should fall off monotonically and rapidly with x. 
So, if one goes to sufficiently large x, the sea quark contribution will be negligible relative 
to the CSV term. One then has to extract the small, x-dependent term in Eq. (|l^) from 
the large term independent of x. As a general rule, the larger the values of x and z at which 
data can be taken, the larger the CSV term will be relative to the 2r-dependent term. 

Note that for a nuclear target, the quark distributions which appear in Eq. ([I^ are 
the nuclear quark distributions, and not the free nucleon distributions. However, for the 
nonstrange quarks a common assumption is that the quark distributions for flavor i, in 
a nucleus with A nucleons, are related to the free ones by qf'^x) = e{x)qf{x). If e(x) is 
independent of quark flavor (most models of nuclear effects on quark distributions make 
this assumption), then this quantity cancels in the ratio in Eq. (|16|), which would then be 
identical to the ratio for free nucleons. 
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III. PREDICTIONS OF CHARGE SYMMETRY VIOLATION IN PION 

LEPTOPRODUCTION 

A. A Simple Model for the Fragmentation Function into Charged Pions 

We want an estimate of the fragmentation function for a quark into a charged pion. Our 
main interest will be to use this simple model to estimate the magnitude of charge symmetry 
violation in this fragmentation function, as this quantity enters into the ratio {z) of Eq. 
(p!7|). Our calculation is based on a method used recently by Nzar and Hoodbhoy pT| . 



They calculated the fragmentation function for a quark into a nucleon plus diquark. The 
fragmentation function can be expressed in terms of the quark field operators and the light 
cone momentum fraction z = P^/k~^ 

4^^P+ = E j ^e-^/^Tr(0|7+^(0)|P;X) (P; X|^(An)|o) (19) 

In Eq. ([TP|), D'^^. {z) is the fragmentation function for a quark of flavor i and four-momentum 
fragmenting into a meson of charge j and four-momentum P^, and the unmeasured state 
X, as is shown in Figure 2(a). We also define null vectors and such that = = 0, 
p ■ n = 1, and n ■ A = for any four vector A. 

To calculate the fragmentation function of Eq. ([T9|) , we need to evaluate the contribution 
for every state X and sum over the complete set of states. The simplest approximation is 
that for charged pions the fragmentation function is dominated by a quark which fragments 
into a pion with momentum and a single quark of momentum {k — P)^; this is shown 
schematically in Figure 2(b). We further assume that we can approximate the amplitude 
for this process as 

(P; {k - P)|^|0) ^ u{k - P)$(P)^^ (20) 

In Eq. (0), a quark with current quark mass m fragments into a pion with mass M and a 
quark with constituent mass niq. The matrix $(A;^) is given by $(A;^) = 0(/i;^)75. 
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Inserting Eq. (^) into Eq. (|T9]), summing over the intermediate states and using the 
kinematic conditions gives 



^^"^""^ 4(1 - zf Jo (27r)2 (p _ ^2y 

= nig + [l — z) m — 2mmg{l — z) . (21) 

The kinematic conditions for the fragmenting quark imply that 

= + ^ 22 

z 1 — z 

while for the form factor, 0, we choose the damping factor 

With these assumptions the integral in Eq. (pT]) can be done analytically, with the result 

/3^^ + -^-A^ . (24) 
z 1 — z 

For our calculations we chose A = 0.4 GeV. For the fragmentation function D^'^ (i.e., an 
up quark fragmenting into a vr^ and a down quark), we chose rriu = 5 MeV and = 353 
MeV. The resulting theoretical fragmentation functions are shown as the solid curve in 
Figure 3(a). We have arbitrarily normalized our fragmentation functions to the experimental 
measurements of the EMC collaboration [|1^,0, at large z. As can be seen, for z < 0.4 there 



is a significant deviation between the calculated and measured fragmentation functions. This 
is not surprising as we expect that the truncation of the complete set of states, X, at a single, 
constituent quark should only be reasonable at large z. 

Using these approximations for the fragmentation function, we can estimate the magni- 
tude of charge symmetry violation in the favored fragmentation function, by recalculating 
{z) and taking into account light quark mass differences. For this fragmentation func- 
tion we used = 8 MeV and rriq = 350 MeV. In Figure 3(b) we plot the charge symmetry 
violating fraction for the fragmentation function, 6D{z), given by Eq. (^). The predicted 
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CSV term is of the order of a few percent. SD{z) < as expected, since involves 
creation of a relatively light uu pair compared to the heavier dd pair created in the case 
of the fragmentation function -D^^. Our simple model predicts a rather large CSV term. 
This is because the pion mass is unusually small, and the light quark mass differences are 
small in magnitude but large as a fraction of the current quark masses. We emphasize that 
this calculation should give only a crude estimate of the size of CSV terms in the quark 
fragmentation function. 

B. Calculations of CSV in Pion Leptoproduction 

In the previous section we proposed measuring the quantity R^{x, z), a ratio of charged 
pion leptoproduction cross sections on the deuteron. As seen from Eqs. (|16D and (0), this 
quantity separates into three parts: one piece which depends only on z, a second which 
depends only on x and is proportional to the quark CSV, and a third part depending on the 
contribution from sea quarks. 

In Fig. [4] we plot the ^-dependent contribution to R^{x,z), i.e. Rf{z) of Eqs. (|T6|) 
and (|13). The solid dots are the values for this quantity using the EMC fragmentation 
functions they do not include a CSV contribution from the fragmentation function, 

and therefore depend only on the ratio of favored to unfavored fragmentation functions. The 
open dots include our estimate of the charge symmetry violating term from the fragmentation 
function, taken from Eqs. (pTsD and (^If). The CSV term is roughly 1% of the total z- 
dependent term. The overall z-dependent term is of order unity; it falls off monotonically 
and smoothly with increasing z, until z ^ 0.6, after which it remains constant and may 
increase somewhat. The solid curve shows this quantity as approximated by the Feynman- 
Field parameterization ||2CI|| . 

In Fig. [5] we plot our predictions for the x-dependent terms in R^{x,z), at = 10 
GeV^. The long dashed curve is the contribution from nonstrange sea quarks to z). 
This depends only on x, and is calculated using the CTEQ3M parton distributions from the 
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CTEQ group |22| . The short dashed curve is our prediction for the parton charge symmetry 
violating term, R^gy{x)] this uses the CTEQ3M parton distributions, plus the bag model 
prediction for valence quark CSV from Londergan et al. |T^. The dot-dashed curve is 
the contribution from strange quarks. This last term is proportional to the strange quark 
fragmentation function. As we do not know this, we have removed the z dependence of 
this term by approximating ^s{z)/ (1 + ^{z)) 1. This will substantially overestimate the 
strange quark contribution, but even with this approximation the strange quarks contribute 
only a small fraction of the total. The solid curve is the sum of the three terms. 

From Fig. [5] we see that for x ~ 0.5, the CSV term is as large as the sea quark 
contribution, and with increasing x the CSV term dominates the x dependence of this ratio. 
We predict the maximum CSV contribution will be of order 0.02 — 0.04. In the small-x 
region dominated by the sea quarks, R^^^^x, z) (at constant z) should decrease rapidly with 
increasing x. However, for x > 0.55 the CSV term is predicted to dominate this ratio, and 
at the largest values of x the terms R^gy{x) + i?£^(a;, z) should be dominated by the CSV 
term. 

If we measure the x dependence of R^{x, z) at constant z, then the x dependent contri- 
bution shown in Fig. [5] will sit on a large and constant z dependent term, estimated in Fig. 
[4]. At small x the sea quark term should be rapidly and smoothly falling with x. The value 
of the constant z-dependent background can be obtained by extrapolating the sea quark 
contribution to zero at large x; the CSV contribution would then be the difference between 
this extrapolated value and the measured value at large x. 

From our previous work on parton charge symmetry violation ITHp^ we predict that 
the two CSV terms in Eq. (|16D (i.e., 6d{x) and 6u{x)) will be roughly equal in absolute 
value (each of them should be of the order 1-2% of the average up -|- down valence quark 
distribution), and they should have opposite signs. As a result, in Eq. (|T6[) the two CSV 
contributions would add constructively, and should produce a term whose value is of the 
order of 0.02 — 0.04. Sather's CSV parton distributions fl^ have the same qualitative 
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behavior as ours. 

We predict that the z-dependent term will be much larger then the x-dependent terms, 
as can be seen from comparing Figs. [4] and [5]. Using the EMC measured values for the 
fragmentation functions, as z goes from 0.4 to 0.8 the ^-dependent term in Eq. (|T6|) should 
vary between approximately 1.5 and 1. So the CSV term is expected to be between 1-4% of 
the 2;-dependent term. 

One major result of our calculations is the prediction that the x and z-dependent parts 
of R^{x, z) will separate. Our current calculation has been carried out in the "naive parton 
model;" we have not included things like higher-order contributions to leptoproduction or 
scaling violations. Since we predict that CSV terms will contribute at the few percent level, 
it will be necessary to investigate whether these higher-order contributions to charged pion 
electroproduction are negligible compared with our lowest-order contributions. Another 
question is whether these additional contributions will preserve the separation of variables x 
and z in the quantity R^{x, z) of Eq. (|T6|). These will require more sophisticated calculations, 
which we are presently undertaking. 

In the HERMES experiment at HERA, the goal is to make precision measurements of the 
spin structure functions. For this one must know very accurately the spin dependence of high 
energy electron scattering from deuterium. This requires precise knowledge of the sources 
of systematic error, so the prospect for obtaining very accurate spin-averaged charged pion 
leptoproduction data is excellent. Only data from deuterium targets is required; efficient 
detection of both signs of charged pions is important, but absolute yields are not required 
as overall normalizations cancel out in the ratio of Eq. ([T6|). 

IV. CONCLUSIONS AND FUTURE EXPERIMENTS 

In conclusion, we have shown that one can use charged pion electroproduction on an 
unpolarized deuteron target to measure the charge symmetry violating [CSV] contribution 
in the nucleon's valence parton distributions. These are expected to be of the order of a few 
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percent, so one has to extract a small term in the presence of a large background. Since one 
of the terms depends only on and the other only on x, one expects to exploit this feature 
experimentally. This would constitute the first direct measurement of charge symmetry 
violation in the valence quark distribution functions. 

In defining parton distribution functions, it is routinely assumed that charge symmetry is 
valid. Surprisingly enough, there are very few tests of this assumption. A reasonable estimate 
of charge symmetry violation would be at roughly the 1% level, but it is conceivable that 
charge symmetry violation could be a few times this value and would not have been observed 
to date. An experimental upper limit for charge symmetry in valence quark distributions 
would be very useful. Measurement of a nonzero effect would be extremely interesting. 
In this paper we have explained why the measurement of charged pion leptoproduction 
from an isoscalar target, for example charged pion electroproduction from deuterium, would 
constitute an excellent test of charge symmetry in the valence quark parton distributions. 
We also estimated the charge symmetry violation in the pion fragmentation functions for 
the first time and showed that it would not interfere with the extraction of information 
about the nucleon valence distributions. As measurements of semi-inclusive charged pion 
production on the deuteron are currently in progress in the HERMES experiment at HERA, 
we are hopeful that suitable data might soon be available. 
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FIGURE CAPTIONS 



1. Fragmentation function for quarks into pions, vs. momentum fraction z. (a) Circles: 
favored fragmentation function triangles: unfavored fragmentation (2;), vs. z. 

Data is from the EMC collaboration, Ref. |T^ for open points and Ref. |T^ for solid points, 
(b) Unfavored/favored ratio A (2;) = i^z) j D'"^ i^z) . Data from the EMC collaboration. 
Solid curve: theoretical prediction of Feynman and Field, Ref. . 



2. (a) Fragmentation of quark with momentum k into a pion with momentum P, and any 
remaining state X. (b) Approximation to fragmentation function where the fragmentation 
is assumed to be dominated by a quark-pion-quark vertex. 

3. (a) Estimate of quark-pion favored fragmentation function D'^J^ {^z\ Solid curve: 
model result given by Eq. (^4]), assuming rriu = 5 MeV, niq = 353 MeV, A = 400 MeV; 



experimental points are those of the EMC group, Ref. [|18],|19|. Dashed curve: fragmentation 
function (z), from Eq. (p^, with = 8 MeV, rriq = 350 MeV. (b) Estimate of charge 
symmetry violation in fragmentation function. The quantity 6D{z) of Eq. ([T^) and the 
model fragmentation functions of Fig. (3a). 

4. The quantity (z) of Eq. (P^). Fragmentation functions taken from the EMC experi- 
ment, Ref. [18,19|. Solid dots: does not include charge symmetry violation in fragmentation 
function; open dots: includes estimate of CSV in fragmentation function, from Eqs. (piSD 
and ^ (and shown in Fig. [3]). Solid curve: the quantity R9{z) using the Feynman-Field 



approximation to the fragmentation functions, Ref. |20]. 

5. x-dependent contributions to R^{x,z) of Eq. (|16]). Long-dashed curve: contribution 



from nonstrange sea quarks; dashed curve: R^gy{x) of Eq. ([T7|), the contribution from CSV 
in nucleon parton distributions; dot-dashed curve: estimate of contribution to R^^^{x,z) 
from strange quarks; solid curve: total contribution including both sea quark and CSV 
contributions. Parton distributions (CTEQ3M) are from the CTEQ group, Ref. at 
= 10 GeV^; CSV term is from bag model calculation of Londergan et ai, Ref. [11^ . 
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